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Abstract. We study 27-day variations of the galactic cosmic ray (GCR) intensity for 2005- 
2008 period of the solar cycle #23. We use neutron monitors (NMs) data corrected and 
uncorrected for geomagnetic disturbances. Besides the limited time intervals when the 27- 
day variations are clearly established, always exist some feeble 27-day variations in the GCR 
intensity related to the constantly present weak heliolongitudinal asymmetry in the heliosphere. 
We calculate the amplitudes of the 27-day variation of the GCR intensity based on the NMs 
data corrected and uncorrected for geomagnetic disturbances. We show that these amplitudes 
do not differ for NMs with cut-off rigidities smaller than 4-5 GV comparing with NMs of higher 
cut-off rigidities. Rigidity spectrum of the 27-day variation of the GGR intensity found in the 
uncorrected data is soft while it is hard in the case of the corrected data. For both cases exists 
definite tendency of softening the temporal changes of the 27-day variation’s rigidity spectrum 
in period of 2005 to 2008 approaching the minimum of solar activity. We believe that a study 
of the 27-day variation of the GGR intensity based on the data uncorrected for geomagnetic 
disturbances should be carried out by NMs with cut-off rigidities smaller than 4-5 GV. 


1. Introduction 

Magnetic cut-off rigidities for galactic cosmic ray (GCR) significantly changes depending on 
the level of disturbances of the Earth’s magnetosphere. According to the Dungey mechanism 
[T] relatively high energy particles rush towards the Earth but are diverted around the Earth 
in circular orbits in the equatorial plane, forming a ring current at several earth radii, which 
causes significant geomagnetic field reduction. These reductions in the terrestrial magnetic field 
strength are measured by the Dst index (disturbance storm time index) [2]. 

An influence of short period changes of the Earth’s magnetic field connected with interaction 
of the Earth’s magnetosphere and disturbed vicinity of the interplanetary space (coronal mass 
ejecta, shocks) on GCR intensity and anisotropy was studied based on corrected and uncorrected 
neutron monitors (NMs) data in [3l HI E] and references therein. At the same time, assessment 
of the 27-day variation of the GCR intensity by NMs data corrected for disturbances of Earth’s 
magnetic field was not studied at all. It is worth to underline that a survey of the 27-day 
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Figure 1. Daily changes of DST index, cut-off variations ARc and the corresponding GCR 
intensity for Kiel {Rc = 2.29 GV), Rome {Rc = 6.32 GV), and Halekala {Rc = 13.3 GV) NMs 
corrected and uncorrected for geomagnetic disturbances for year 2005. 


variation of the GGR intensity is fundamental to understand the influence of the heliolongitudinal 
asymmetry of solar wind and solar activity on behavior of the Earth’s magnetosphere and 
atmosphere, and its consequences for the global climate processes. In order to address this issue 
an aim of this paper is: (1) to estimate how the changes in the cut-off rigidities of various NMs 
influence the amplitudes of the 27-day variation of GGR intensity and (2) to study changes in 
rigidity spectrum of the 27-day variation of the GGR intensity in 2005-2008, the time interval 
of the unusually prolonged deep minimum in the solar activity 7^23. 

The recent solar minimum 23/24 provided a unique opportunity to study recurrent variations 
of the GCR intensity under relatively stable conditions. Recurrent variations connected with 
corotating structures (~ 27 days), at the end of 2007 and almost for the whole year 2008 were 
clearly established in all solar wind and interplanetary parameters. Consequently, the 27-day 
recurrent variations of cosmic ray intensity were clearly visible in a variety of cosmic ray counts 
of neutron monitors (e.g., [HE]) and space probes (e.g., [8] ). 

2. Analysis of NMs data corrected and uncorrected for geomagnetic disturbances 

To find changes of cut-off rigidities ARc for different neutron monitors during the analysed 
period of 2005-2008 we have employed the spectrographic global survey (SGS) method devel¬ 
oped by [9] . Using ground-based measurements of GCR from the worldwide network of stations 
the SGS method allow to collect information on energy and pitch-angle distribution of primary 
GCRs in the interplanetary magnetic field (IMF), as well as on variations in the planetary sys¬ 
tem of geomagnetic cut-off rigidities per each observation hour. This method provides also the 
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Figure 2. The same as in Fig. [T]but for year 2006. 


possibility to use the whole set of ground-based recording equipment (the worldwide network 
of neutron monitors at any level of the Earths atmosphere, ground-based and underground me¬ 
son telescopes, etc.) for analysis. Along with phases of the first and second harmonics of the 
pitch-angle anisotropy, it can determine the rigidity spectrum of the isotropic component and 
anisotropy, obtain information about the IMF orientation from data on the phase of pitch-angle 
anisotropy. Moreover, the SGS method can determine variations in the GCR planetary system 
per each observation hour or per shorter time intervals when the CR intensity near Earth goes 
up during solar proton events. 

The essence of the SGS method is as follows: if amplitudes of variations in primary GGR 
intensity AJ/J are small, then variations in intensity of secondary particles AI/I are connected 
with AJ/J by the following relation: 

= - / ARc{a,P)W/[{Rc{a,/3),a,Xhi]sinj3d/3da+ ( 1 ) 

-^c X Jo 

/•Stt roo A 7 

+ — —{R,a, fJihi)sinj5dj3dadR 

Jo Jo JRc(oi,/3) J 

ji 

where, -jr{hi) is the amplitude of variations in flux of secondary particles of type i observed 
at geographical point c at height hi in the Earth’s atmosphere at time t; ^{R,a,fJ) 
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Year 2007 


Figure 3. The same as in Fig. [T]but for year 2007. Due to shutdown of the Haleakala station 
the analogous changes for the Potchefstroom NM (i?c = 7.10 GV) are presented. 


are the variations in particle intensity at the boundary of the atmosphere at the given 
point; a, (3 are the azimuth and zenith angles of primary particles’ arrival to the boundary 
of the atmosphere; R is the magnetic rigidity of particles; Rc{a,j3) is cut-off rigidity in 
directions (a, /?), and ARc{a,/3) are its possible variations. The W^{R,a,l3,hi) is the 
function determining a connection between ^^{hi) and ^{R, a, j3) satisfying the normalisation 
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condition: Jq^ da /g^ sin(3dp p^W^{R,a, /3,hi)dR = 1. Assuming that direction of the 
particle motion in the given coordinate system is determined by the angles and A, and the 
IMF direction by the angles ipo and Aq, then the particle pitch angle 9 can be expressed by 
following characteristics /r = cos9 = sinXsinXo + cosXcosXocosltp — ipo). Simultaneously, we can 
represent the function A) as a series V’, A) = J2^=o^n{R)PnifJ') where Pn{lj) are 

the Legendre polynomials. Taking into account the simplifications and restrictions presented in 
detail in m we obtain 

A ji 10 Tn-n 

= -AR,Wi{Rc, + E E (2) 

C j=0k=l 

where Ag^ = aok', ^ik = oifcCos'i/’ocosAo, A 2 k = aiksin'ipocosXo] A^k = aiksinX^, 
A^k = a2kCos'^'ijjocos'^Xo, A^k = a2kCos'ilJosimlJocos'^Xo, A^k = a2ksin‘^i’ocos'^Xo, A-jk = 
a 2 kCos'il)QCOsXQsinXQ, Agk = a 2 ksimpQCOsXosinXo, Agk = a 2 ksi'n?Xo, AiQk = a 2 k- The coefficients 
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Figure 4. The same as in Fig. [3] but for year 2008. 


are determined based on the asymptotic angles derived from calculations of the particles 
trajectory m- Summarizing, the problem reduces to estimation of the following unknown 
parameters: Ai?c-showing the variation of cut-off rigidity; aoi, ao 2 , aosr oii) ai 2 , 021 , ^22 -showing 
amplitudes and rigidity dependencies of the zero, first and second harmonics of the pitch-angle 
GCR distribution; 'i/jq and Aq - showing the IMF direction averaged over the Larmor circle 
(accurate to vr). To solve the problem in hand, it is essential that the network of stations and 
detector arrays at stations used in the calculations provide abundance of the system of above 
mentioned equations and linear independence of coefficients There should be at least two 

detectors with different Wl{Rc, hi) at the stations where GCR variations are determinable. The 
details of the SGS method are presented in m- 

Applying the SGS method we have calculated the daily changes of the ARc by means of 27 
neutron monitor stations hourly data during years 2005-2006 versus the cut-off rigidity Rc on 
10 Jun 2004, and for period of 2007-2008 versus Rc on 5 November 2007. Figs. [T]- [J] present 
daily changes of Dst index, cut-off variations ARc and the corresponding GCR intensity for 
Kiel, Rome, Potchefstroom and Haleakala NMs corrected and uncorrected for geomagnetic 
disturbances for 2005 (Fig. [T|), 2006 (Fig. [2|), 2007 (Fig. [3|) and 2008 (Fig. [3|). Fig. HJ-Fig- El 
show that the corrected and uncorrected GCR intensity do not differ much for NMs with cut-off 
rigidities smaller than 4-5 GV, while for NMs with higher cut-off rigidities the difference is as 
notable as cut-off rigidity increases i.e. a level of changes of cut-off rigidities is larger. One can 
see that the correction for geomagnetic disturbances is the largest when the rapid changes of Dst 
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Figure 5. Daily GCR intensity corrected and uncorrected for geomagnetic disturbances with 
filtered 27-day wave and the corresponding correction for Kiel and Rome NMs in 2007-2008. 


index are observed (related with shocks and Forbush decreases), but becomes weaker in quite 
stable conditions. 

In this paper we consider in detail 27-day wave in corrected and uncorrected GCR 
intensity and the effect of geomagnetic correction on periodic variations. Daily GCR intensity 
corrected and uncorrected for geomagnetic disturbances with filtered [12] 27-day wave and the 
corresponding correction for Kiel and Rome NMs in 2007-2008 are presented in Fig. [SJ As an 
example Fig. [S] presents only few Bartel’s rotations in 2007-2008 with excellent quasi-periodic 
changes (~ 27 days) in the GCR intensity. 

3. Rigidity spectrum of the 27-day variation based on the corrected and 
nncorrected NMs data 

Usually, a power law rigidity R spectrum of any classes of the GCR intensity variations is 
characterized by an exponent 7 , as ~ R~'^. The exponent 7 is an important instrument to 
study a dependence of amplitudes of the GCR intensity changes on energy of primary GCR 
particles. The value of 7 is determined by the structure of the IMF turbulence responsible for 
the scattering of GCR particles to which neutron monitor respond (e.g. [laiiaiie]). A character 
of this dependencies is different for various classes of the GCR intensity variations. In [18] was 
shown that temporal changes of the exponent 7 for the 11-year variations |13| inversely correlates 
with the changes of the exponent 7 for the 27-day variations of the GCR intensity. 

In this paper the values of the rigidity spectrum {6D{R)/D{R)) exponent 7 of the 27- 
day variation of the GCR intensity were calculated using both, corrected and uncorrected for 
geomagnetic disturbances data of NMs with different cut-off rigidities. 

We have calculated the amplitudes of the 27-day variation of the GCR intensity by means of 
daily data using the harmonic analyses method e.g. m during each Bartel’s rotation period: 

N/2 

x{kAt) = — ''^[anCos{—nkAt) + bnSin{-^nkAt)], 

^ n=l ^ ^ 


( 3 ) 





































































Date, Years 2005-2008 


Figure 6. (Top panel)The temporal changes in the amplitudes of the 27-day variation of the 
corrected and uncorrected GCR intensity for Kiel and Rome NMs; (Bottom panel) The exponent 
7 of the rigidity spectrum of the 27-day variation of the corrected and uncorrected GCR intensity 
in 2005-2008. 


where x{kAt) designates the daily data of the GCR intensity, N = 27 and the coefficients have 
a form: oq = ^J2k=ix{kAt), an = jfJ2k=i[x{kAt)cos^^], bn = jfJ2k=i[x{kAt)sin^^]. 
Applying the above presented formulas we determine the amplitude of the 27-day variation 
(A27I) of the GCR intensity as A27I = a\ + h\. 

The rigidity spectrum of the quasi-periodic variation was calculated by means of the 
amplitudes of the 27-day variation of the GCR intensity based on the method presented, e.g. in 

mm- 


SDjR) ^ I AR-^ R < Rmax .4^ 

D{R) [0 R> Rmax ^ ’ 

Rmax designates the upper limiting rigidity beyond which the quasi-periodic variation of the 
galactic cosmic ray intensity disappears (100 GV). Detailed description of the 7 calculation is 
described in [19]. 

An investigation of the long-period changes of the power law rigidity spectrum of the 27-day 
variation of the GCR intensity (e.g. in the period 2005-2008) requires a data from neutron 
monitors with a long term stability. Unfortunately, only a few neutron monitors from worldwide 
network of stations satisfy this requirement. For this reason, we chose Kiel and Rome neutron 
monitors data to calculate the rigidity spectrum of the amplitudes of the 27-day variation of the 
GCR in the period 2005-2008. To use two neutron monitors for the calculation of the reasonably 
reliable rigidity spectrum exponent 7 there must be considerable difference between the cut-off 
magnetic rigidities {Re)- This demand is satisfactorily fulfilled for Kiel and Rome neutron 
monitors; for Kiel neutron monitor R^ = 2.29 GV and for Rome neutron monitor Rc = 6.32 GV 
(or the corresponding median rigidity of response for Kiel neutron monitor equals 17 GV and 
for Rome - 23 GV, respectively. A validity of this approach was shown in papers usmsi. 

Fig. [ 6 ] (top panel) demonstrates the temporal changes in the amplitudes of the 27-day 















































variation of the corrected and uncorrected GCR intensity for Kiel and Rome NMs; and (bottom 
panel) the matching exponent 7 of the rigidity spectrum of the 27-day variation of the corrected 
and uncorrected GCR intensity in 2005-2008. Unfortunately, in the first half of the year 2006 
calculation of the exponent 7 based on the corrected data was not possible due to lack of the 
discrepancy between the amplitudes of the 27-day variation for the used NMs. Fig. [ 6 ] shows that 
rigidity spectrum of the 27-day variation of the GCR intensity found by uncorrected data of 
NMs is soft while it is hard for corrected data. For both cases there is visible a clear tendency of 
the softening of the temporal changes of the rigidity spectrum from 2005 to 2008 (approaching 
minimum of solar activity) coinciding with results obtained in |16j . 


4. Summary 

(i) Comparison of the amplitudes of the 27-day variation of the GCR intensity calculated based 
on the NMs data corrected and uncorrected for geomagnetic disturbances showed that these 
amplitudes do not differ for NMs with cut-off rigidities smaller than 4-5 GV in comparison 
with NMs of higher cut-off rigidities. This indicates that to study the 27-day variations of 
the GCR intensity by uncorrected data of NMs it is necessary to use data of NMs with cut 
of rigidities smaller than 4-5 GV. 

(ii) Rigidity spectrum of the 27-day variation of the GCR intensity found by uncorrected data 
of NMs is soft, while it is hard for data corrected on geomagnetic disturbances. For both 
cases, there is an apparent tendency of softening temporal changes of the rigidity spectrum 
from 2005 to 2008 period approaching minimum of solar activity. 
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